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We present an analytical model which permits the calculation of effective material parameters 
for planar metamaterials consisting of arbitrary unit cells (metaatoms) formed by a set of straight 
wire sections of potentially different shape. The model takes advantage of resonant electric dipole 
oscillations in the wires and their mutual coupling. The pertinent form of the metaatom determines 
the actual coupling features. This procedure represents a kind of building block model for quite 
different metaatoms. Based on the parameters describing the individual dipole oscillations and 
their mutual coupling the entire effective metamaterial tensor can be determined. By knowing these 
parameters for a certain metaatom it can be systematically modified to create the desired features. 
Performing such modifications effective material properties as well as the far field intensities remain 
predictable. As an example the model is applied to reveal the occurrence of optical activity if the 
split ring resonator metaatom is modified to L- or S-shaped metaatoms. 
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I. INTRODUCTION 

Metamaterials extend the optical effective material re- 
sponse of natural media. They are artificial materials 
that allow to tailor light propagation properties by a care- 
ful design of the mesoscopic metaatoms the metamaterial 
is made of. By controlling the geometry and selecting the 
material dispersion of the metaatom, novel effects such 
as negative refractioiP"^, optical cloaking^ ^, as well as 
a series of optical analogues to phenomena known from 
different disciplines in physics have been observed ^"^H^. 

In principle, the material response is mathematically 
best described in terms of constitutive relations; lead- 
ing to tensors that relate the electric displacement and 
the magnetic induction to the electric and magnetic 
fields, respectively. In the initial stage of research on 
metamaterials emphasize was put on exploring materi- 
als that potentially lead to a bi-axial anis otrop ic (lin- 
ear dichroism) effective material response^^^^^^^'^. Re- 
cently research was also extended towards the explo- 
ration of metaatoms that affect off-diagonal elements of 
the effective material tensors (elliptical dichroism). It 
expands the number of observable optical phenomena, 
leading to, e.g., optical activity"^^, bidirectional and 
asymmetric transmission^^ or chirality-induced nega- 
tive refractiorP^H^. 

In our understanding optical activity comprises all ef- 
fects on light propagation resulting from the nonlinearity 
of the polarization eigenstates, hence including phenom- 
ena like circular dichroism, all effects of genuine three- 
dimensional chirality and the optical manifestations of 
planar chiralit}'^^. 

However, despite the option to resort to rigorous com- 
putations for describing the light propagation on the 



mesoscopic level of the metaatoms, an enduring issue in 
metamaterial research is the question how the effective 
material tensor looks like for a certain metamaterial. 

In general, investigating the geometry of the metama- 
terial (the metaatoms geometry and their arrangement) 
allows to determine the form of the effective material 
tensors in the quasi-static limit as extensively discussed 
in Ref. [191 From such considerations it is possible 
to conclude on the symmetry of the plasmonic eigen- 
modes sustained by the metaatoms and on the polar- 
ization of the eigenmodes allowed to propagate in the 
effective medium^^. But in order to determine the actual 
frequency dependence of the tensor elements, more ex- 
tended models are needed which start in their description 
of the metaatom properties from scratch^^. Such models 
are required to be universal, simple and assumption-free 
to the largest possible extent. 

Here we outline an approach which meets these re- 
quirements. It is based on conceptually decompos- 
ing the complex metatoms into a set of coupled plas- 
monic entities that sustain the excitation of dipolar 
resonances^^. The knowledge of the plasmonic proper- 
ties of these dipoles and their coupling suffices to derive 
the material response and the symmetry of the eigen- 
modes. This, in turn, permits to predict the observable 
quantities in the far-field, such as the polarization and 
frequency-dependent reflected and transmitted complex 
amplitudes. 

The most appealing aspect of our model is that once 
the plasmonic entities and their coupling strengths are 
characterized, far-field properties remain predictable by 
the analytical model even if substantial modifications of 
the metaatom geometry have been made. Even a mod- 
ification that leads to a different symmetry of the ma- 
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terial tensor does not prevent a quantitative description 
of the dispersive behavior of the tensor elements. To be- 
come specific, we start our investigations with an optical 
inactive, bi-axial anisotropic metamaterial, namely the 
split-ring resonator (SRR). From the observable far-field 
quantities we derive the properties of the dipolar oscil- 
lators that govern the plasmonic response of the three 
wires forming the SRR, namely the individual eigenfre- 
quencies, oscillator strengths and damping constants as 
well as their mutual coupling strengths. By relying on 
these quantities, we determine the effective properties of 
metamaterials consisting of modified metaatoms with re- 
spect to the initial SRR. 

We focus on two modifications that evoke optical ac- 
tivity. With our analytical model we predict the effec- 
tive properties of these metamaterials. We then compute 
the optical coefficients for a slab made of these metama- 
terials and compare it to rigorous simulations. Excel- 
lent agreement between the analytically predicted and 
the rigorously calculated optical coefficients is observed 
throughout our work. Therefore, the proposed method 
can be used for the parameter retrieval without resorting 
to rigorous simulations. And since it is based on analyti- 
cal considerations the approach potentially allows a large 
variety of metaatom modifications and to systematically 
tailor its effective properties. Hence our approach pro- 
vides a powerful and versatile tool for a systematic anal- 
ysis of achievable material properties by varying only a 
few constituents that may couple in some well-defined 
ways. 

Furthermore, we will also show that such an analyti- 
cal treatment provides further insight into metamaterial 
properties. Specifically, we show that it is possible to 
directly infer that the model's predictions are valid in 
terms of the Casimir-Onsager relations^^ the require- 
ment for time-reversal and reciprocity in linear medisP^. 

Thus, the main aspects of the paper can be summa- 
rized as follows. At first, the localized carrier dynamics 
occurring in metaatoms may be properly described by 
a set of coupled oscillators, representing the decomposi- 
tion of the metaatom in nanowire pieces. Second, the 
dynamics of these oscillators, determined by the shape of 
the nanowires and their coupling, result in electric dipole 
polarizabilities that permit the calculation of the effec- 
tive permittivity tensor. The main advantage of this 
simple approach is that modifications of a metaatom, 
for which the oscillator dynamics and parameters have 
been found, leave the effective permittivity tensors pre- 
dictable. Hence, the far field refiectance and transmit- 
tance can be calculated. We will show that this holds 
also for modifications changing the character of the eigen- 
states from linear to non-linear polarized. Moreover, the 
approach can be useful to determine the effective per- 
mittivity tensor for optically active metamaterials, whose 
eigenstates are in general elliptically polarized, by inten- 
sity measurements of linearly polarized light. 
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FIG. 1. (Color online) (a) The original SRR structure (left), 
the first modification, namely the L- (center) and the second 
modification, the S-structure (right), (b) The SRR together 
with the carrier oscillators, marked by black dots, which are 
used to phenomenologically replace the SRR. 

II. THE PLANAR SRR METAATOM 

To reveal the versatile character of our approach we 
will directly start by conceptually replacing the planar 
SRR geometry, shown in Fig. [ijb), by a set of coupled 
oscillators. Each of the oscillators introduced here repre- 
sents a metaatoms piece, i.e. a straight nanowire that is 
coupled to its electrically conductive neighbors. The os- 
cillators accounting for the isolated nanowire are associ- 
ated with the excitation of carriers representing the free- 
electron gas of the metal. Hereby the carrier dynamics 
is solely influenced by the external electromagnetic fleld, 
including contributions from adjacent nanowires. The 
dynamics becomes resonant for the eigenfrequency of the 
localized plasmon polariton resonance of the individual 
nanowire. This localized resonance described by the indi- 
vidual oscillator corresponds to the fundamental electric 
dipole mode since the dimensions of the nanowires form- 
ing the metaatom investigated here are small compared 
with the wavelength. We will show later that the carrier 
oscillations of metaatoms assembled from several coupled 
wires still correspond to electric dipole polarizabilities. 
This holds as long as the wires are assembled in-plane. 
For out-of-plane structures higher order multipoles come 
generally into pla}'^^ 

Here the spatial coordinate represents the elongation of 
a negatively charged carrier density driven by an external 
electromagnetic fleld, which is the usual assumption in 
plasmonic^^. 

We show below that this assumption is sufficient to 
entirely predict the optical response. The equations for 
the coupled oscillators are 

9^ do qi ^ 

g^xi + 71— + ^01^1 + ^21^2 = --^x^ 

+ 73 ^^3 + OJq^Xs - (J23y2 = ""^^^ 

(1) 

where we have assumed nearest neighbor coupling be- 
tween adjacent, i.e. conductively coupled dipolar oscil- 
lators. The oscillators are driven by the electric fleld 
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component of an external illuminating field propagating 
in z direction (normal incidence is assumed throughout 
the entire manuscript). Excitation of the oscillat ors by 
magnetic field components can be safely neglectecP^'^. 
In Eqs. ^ 

is the damping, uoj the eigenfrequency, aij the cou- 
pling constant and qj is the charge for the three oscilla- 
tors (i^j) G (1,2,3), respectively, similarly to those in- 
troduced in Ref. 29, The coordinates (xi, 7/2,^3) them- 
selves are understood as the displacement of the nega- 
tively charged carriers representing oscillating currents or 
multipole moments where both approaches can be used 
to determine the effective material response of an artifi- 
cial medium exhibiting a carrier dynamics as described 
by Eq. Here we apply the multipole approach to 

map the displacement onto electric dipole moments, since 
this is consistent with our previous work4^^^. 

Accounting for the contributions of the electric dipole 
and quadrupole as well as the magnetic dipole moment 
the wave equation reads as 

AE(r,a;) = — /ioeo^^E(r, cj) — uj'^ j^qP {r ^ uj) 

+cc;^/ioV • Q(r, u) - iu/j^oV x M(r, cj), (2) 

We put emphasis on the fact that the higher order mul- 
tipole moments, leading to Q,M, appear in general^'', 
but do not provide any contribution in the present con- 
figuration, as mentioned above. Thus, there is no ef- 
fective magnetic response (effective permeability tensor) 
as expected and it is sufficient to consider the disper- 
sion in the susceptibility, resulting in a linear effective 
permittivity tensor onl}^^^. This will hold for all pla- 
nar configurations with an illuminating field invariant in 
the X — 7/-plane, hence no spatial dispersion occurs which 
would result in an artificial magnetic response. Substi- 
tuting the displacements ([T]) into the definition of the 
dielectric polarization^^ we can introduce the effective 
susceptibility tensor x(^) 

N 



1=1 



(3) 



where r] accounts for the carrier density. This effective 
susceptibility can be easily calculated. As usual we define 
the effective permittivity tensor as 



■x(w), 



(4) 



that governs the wave propagation of an incident 
plane wave in an effective medium composed of SRR 
metaatoms. 

Next we consider the possible eigenmodes of Eqs. ([T]) 
for the two polarization directions {x or y). It suffices 
to investigate these two polarizations as long as the sys- 
tem under consideration is linear. In order to describe 
a SRR with two identical side arms we set cjqi = oJos = 

^Ox, 7l = 73 = Ix, Cr21 = Cr23 = cr and Ql = Qs = -Qx- 

For the oscillator associated with the SRR base having 



a different geometry, we put 6^02 = ^Oy^ 72 = "fy and 
q2 = —Qy. We note that the latter distinction could have 
been dropped if the geometry of all constituents of the 
SRR would have been the same. For the reasons obvious 
from the consideration below we refrain from doing so. 

For a polarization of the incoming plane wave paral- 
lel to the X-axis we can solve Eqs. ([T]) and obtain the 
following displacements in Fourier domain: 



y2{uj) = 0, 



(5) 



where we have introduced Ax{uj) = — cj^ — iujx- 
For the polarization in ^/-direction we obtain 



xi{uj) = -xs{uj) 



m Ax{uj)Ay{uj) — 2cr^ 



■Ey{Z,Uj), 



^ m Ax(uj)Ay{uj) -2(j^ ' ^' 



(6) 



with Ay(uj) = ujQy — uj^ — iuj^y. 

Considering the eigenmodes for x-polarization, we ob- 
serve that two parallel dipoles {xi = xs) are induced, 
while due to symmetry constraints no dipole is induced 
in ^/-direction {y2 = 0), see Fig. [2|c). By contrast, be- 
sides a dipole in ^-direction the ^/-polarized illumination 
induces oscillating dipoles in x-direction in both arms. 
But due to the anti-symmetric oscillation X3 = —Xi [Eq. 
(|6|], the dipoles in the SRR arms [Fig. |2|d)] do not radi- 
ate into the far-field because they oscillate tt out-of-phase 
and interfere destructively. Hence, no cross-polarization 
is observed and the far-field polarization equals that of 
the illumination. Later we will prove that any radiation 
emerging from cross-polarized dipole moments will result 
in an optical activity, as expected. 

By substituting Eqs. ([sj |6| in Eqs. (|3| we get the 
susceptibility tensor for the pertinent SRR configuration 



x(^) = 



Xxx 






XyyH 





meo Aj;{uj)' 
^^^^"^^ meoA,{uj)Ay{uj)-2a^' 



(7) 



with the polarization P = — 7^ {Qxi^i + X3), g^?/2, 0)^, 
according to Eq. As expected, the susceptibility 

tensor is diagonal. Hence, the eigenmodes of the ef- 
fectively homogenous medium are linearly polarized and 
orthogonal. Due to the polarization dependent carrier 
dynamics the SRR shows a linear dichroitic behavior. 
Hence, our model correctly predicts the linear polariza- 
tion eigenstates as required by the mirror-symmetry with 
respect to the x — z-plane. The only unknown parame- 
ters are those describing the oscillators and their cou- 
pling strengths. They can be determined by matching 
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FIG. 2. (Color online) The rigorously calculated (FMM) far- 
field transmission/reflection spectra compared with those ob- 
tained by the analytical coupled dipole model (DM) for the 
two indicated polarization directions [(a) x-polarization and 
(b) ^/-polarization] . The stationary carrier elongation (nor- 
malized imaginary part of xi,3 and ^2) for the two correspond- 
ing polarizations [(c) and (d)]. (e),(f) The exactly retrieved 
parameters compared to the analytical calculations. 



the optical coefficients of an effective medium whose per- 
mittivity is described by Eqs. ^ to the spectra obtained 
by rigorous simulations or far-field measurements of the 
structure. 

In order to find the oscillator parameters, we 
performed numerical Fourier Model Method (FMM) 
calculations^^ of a periodic array of gold SRR's^ similar 
to those reported in Ref. 2, These numerical far-field 
observables (reflected and transmitted intensities) were 
fitted by using the effective permittivity tensor Eq. Q 
for both polarization directions in a conventional transfer 
matrix formalism that computes reflection and transmis- 
sion from a slab of equal thickness^^. Results are shown 
in Fig. [ga) and [lib). 

Once the unknown parameters have been found the 
frequency dependent stationary elongations of the oscil- 
lators can be determined as shown in Figs. |2|c) and 
(d). They can be used to identify the carrier oscillations 
of the different plasmonic eigenmodes. These carrier os- 
cillations are shown on top of Fig. ([2| and their hor- 
izontal position relates to the respective resonance fre- 
quency. The observed eigenmodes are documented in 
literature^^. Figures [2|e) and (f) show the rigorously 



retrieved effective material properties together with the 
ones of the analytical model. Excellent agreement is ob- 
served. The rigorous results were obtained by applying 
a common parameter retrieval based on an inversion of 
the matrix formalism to calculate the effective parame- 
ters on the basis of complex reflected and transmitted 
amplitudes for a certain slab thickness. We underline 
that the unknown oscillator parameters can be obtained 
by comparing the far-field intensity only. The derived 
effective material parameters are in excellent agreement 
with the rigorous results, but they were derived without 
the necessity of knowing the complex- valued fields. This 
will be advantageous for optical active structures where 
the experimental determination of these parameters is in 
general involved (e.g. phase resolved measurements are 
not required). 



III. THE L - METAATOM 

As outlined in the introduction low-symmetry 
metaatoms are investigated in order to reveal new op- 
tical phenomena like asymmetric transmission. We will 
therefore extend the previous considerations towards op- 
tically active effective media by rearranging the SRR con- 
stituents. In the following we shall rely on the oscilla- 
tor parameters obtained by the fitting procedure above. 
Using these parameters in conjunction with the analyti- 
cal expressions for the rearranged constituents as derived 
below, the optical response of the modified metaatoms 
can be predic ted. T he first investigated structure is the 



L - metaatonj ^^ l ^^ l In this geometry one of the SRR 
arms is omitted in order to prevent the cancelation of 
the far-field that originated from the antiparallel dipole 
moments. Hence, it is expected to obtain polarization 
rotation. To get specific in Eq. ([T]) one horizontal oscil- 
lating dipole, e.g., dipole '3', has to be dropped. Hence, 
we obtain for x-polarization 



y2 



m Ax{uj)Ay{uj) — cr^ 



m Ax{uj)Ay[uj) — a'^ 
and for ^/-polarization 



(8) 



V 

y2 



m Ax{uj)Ay{uj) — a' 
qy Ax{uj) 



;Ey{z,Uj), 



m Ax{uj)Ay{uj) — cr^ 
From the resulting polarization Eq. ([3| 



Ey{Z,Uj). 



J 



P^{z,uj) = -r] I qyyl \ , j e [x,y] 



(9) 



(10) 
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FIG. 3. (Color online) The far- field response of the L - struc- 
ture for X — (a) and y — polarization (b) . In addition to 
the numerical (FMM, spheres) and the fitted data (DM, solid 
lines) the predicted spectra incorporating the SRR parame- 
ters (dashed-dotted lines) are plotted. (c),(d) In contrast to 
the SRR both eigenmodes can be excited for each polarization 
direction. The respective numerical cross-polarization contri- 
butions (FMM, circles) compared with the model predicted 
(dashed-dotted lines) and the fitted (solid lines) values are 
shown in (e) and (f). Note that figures (e) and (f) are identi- 
cal as required for such kind of effective media and are only 
shown for completeness. 



the susceptibility reads as 



Xxx 

Xyx{oj) Xyy{^) 





Xxx{(^) 

Xyy{(^) 
Xyx{^) = Xxy{^) 



Ay{uj) 



mco Ax{uj)Ay{Lj) 



meo Ax{uj)Ay{uj) - 

QxQyr] cr 

meo Ax{uj)Ay{uj) - c 



(11) 



The most significant change compared to the SRR is 
the appearance of off-diagonal elements in the suscepti- 
bility tensor x{uj) which is, however, symmetric leading 
also to Sij{uj) = eji{uj). This symmetry relation is impor- 
tant because it is required for time reversal, known as the 
Onsager-Casimir principl^^^''^. As expected the tensor 
of the effective permittivity has t he same form as that for 
planar optical active medis^^ resulting in asymmetric 



transmission due to elliptical dichroism. Note that the 
optical response would change dramatically if both arms 
would be identical. In this case the diagonal elements 
Xii{^) are identical too and the tensor can be diagonal- 
ized by an rotation of 7r/4. So the polarization eigenstates 
would be linear and the effective medium would be lin- 
early dichroitic. This is clear since the metaatom would 
have an additional mirror symmetry with respect to the 
plane defined by the surface normal and the line x = 
see Ref. \M 

Another difference while comparing the L - to the SRR 
metaatom is that the splitting a between both resonances 
is reduced (by a factor of >/2), which follows from drop- 
ping one of the SRR arms. 

In order to check whether our simple description is 
valid and to reveal the relation between the SRR and the 
L - structure eigenmodes, we performed again numerical 
FMM simulations for the L - metaatom and compared 
the results to the analytical model in Fig. [S] Since 
we deal now with more complex effective media where 
the full tensorial nature of the permittivity tensor has to 
be taken into account the standard transfer matrix algo- 
rithm cannot be applied. Hence, it is challenging to de- 
termine the scattering coefficients analytically' or even 
to invert them to retrieve effective parameters directly. 
Therefore we used an adapted Fourier Modal Method to 
determine the transmitted and reflected intensities^^! In 
a first approximation we used the parameters which we 
determined for the SRR. Based on these parameters we 
calculated Xij{^) the far-field intensities as shown 
in Fig. [3]^a,b) for the L - structure. The associated 
eigenmodes for the carriers are shown in Fig. |3|c,d). 
The curves for both the co-polarized Fig. |3|a,b) and the 
cross-polarized intensities Fig. [3]^e,f) are in good qual- 
itative agreement for the rigorous FMM results (dotted 
line) and the analytical model (dashed line) based on 
the previously derived SRR oscillator parameters. Al- 
though there are deviations between the actual resonance 
strength, the agreement, e.g. for the resonance positions, 
is obvious. Note that the intensities, in particular for the 
cross-polarized fields, can be solely predicted by the coef- 
ficients obtained from the SRR. There, an almost perfect 
agreement is observed. 

In a second step we adapted the oscillator parame- 
ters in order to fit the exact calculations [solid lines in 
Fig. [3]^a,b,e,f)] providing an almost perfect coincidence 
with the numerical values^^. Note that the fitting is done 
only for the co-polarized intensities from which the cross- 
polarized intensities follow. 

A last step yields the effective permittivity tensor, that 
is inherently accessible and already applied in order to fit 
the spectra in Fig. |4] It can be seen that for the two po- 
larization directions two eigenmodes appear as Lorentz- 
shaped resonances for the effective permittivities. They 
differ in strength, due to the different geometrical pa- 
rameters of the L - arms. Fig. |4ja,b). The off-diagonal 
elements €ij{uj) are identical as discussed before, [Fig. 
|4jc)]. Considering especially the second resonance for 
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FIG. 4. (Color online) The diagonal effective permittivity tensor elements of the L - metaatom: (a) e(c<;)xx and (b) e(c<;)yy. The 
arrows indicate the current flow, given by direction of the carrier oscillation for the particular eigenmode. (c) The off-diagonal 
elements e{(jj)ij comprising a Lorentz resonance around i> — 7,000 cm~^ and an anti-Lorentz resonance at P = 11,000 cm~^ . 



the off-diagonal tensor elements near v = 11,000 cm~^ , 
we observe a Lorentzian anti-resonance that might sug- 
gest gain within the system due to the negative imaginary 
part. However, the nature of this resonance can be ex- 
plained by the introduced formalism as well. Since the 
permittivity is proportional to the susceptibility Eg. Q 
and hence also to the carrier displacements Eqs. (11), 
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a negative sign corresponds to a phase difference of tt 
between both oscillating carrier densities (antiparrallel 
oscillations), while for the positive Lorentz resonance at 
around V = 7,000 cm~^, both are oscillating in-phase 
[sketched by the arrows in Fig. ga)]. Thus, the nega- 
tive sign in the imaginary part of the permittivity can by 
fully explained by means of the mutual interplay of the 
coupled oscillators. 



IV. THE S - METAATOM 

In order to verify the model we are going to inves- 
tigate another modification of the SRR, namely the S 
- structurd^^^ [Fig. [ija)]. To observe optical activity 
with the same number of coupled entities as for the SRR 
its mirror symmetry has to be broken. Therefore one of 
the SRR arms (e.g. xs) is turned with respect to the 
SRR base. 

This opening of the SRR structure is expected to en- 
able the observation of two modes for x-polarization, 
since the oscillator in the base (^2) can now oscillate in- 
phase or out-of-phase with the two remaining x-oriented 
oscillators of the S - structure that are excited by the x- 
polarized electric field. For ^/-polarization the situation 
is similar, but now the two oscillators in the horizontal 
arms (xi, xs) are allowed to oscillated in-phase or out- 
of-phase to the excited oscillator in the base {1/2). Math- 
ematically this modification can be considered by setting 
(723 = —o'2i = CF (see Eq. [T]), while all other parameters 
appear similar to the ones applied for the SRR. 

With these initial assumptions, which reflect all modifi- 
cations to the geometry, the calculations can be repeated 
in analogy to those for the SRR and the L - structure. 
Thus, we obtain the elongations for x-polarized excita- 



m 



m Ax{uj)Ay{uj) — 2cr^ 
and for ^/-polarized excitation 



(12) 



as well as the respective effective susceptibility tensor 



m . 



= 

Xyy{^) = 
Xyx{^) Xxy{^) = 



Xxx{^) Xxy{^) 
Xyxi^) Xyyi^) 





2Ay{u) 



meo Ax{uj)Ay{uj) — 2<j^ ' 

meo Ax{^)Ay{^) — 2cr^ ' 

qxQyV 2a 

meo Ax{^)Ay{^) — 2cr^ 



• (14) 



In Eqs. (14) we have used the polarization Eq. ([3| 
which is found to coincide with that of the SRR struc- 
ture, whereas the elongations are different. With respect 
to the splitting of the resonances, we expect the same res- 
onance positions as for the SRR for ^-polarization, since 
2a appears in the denominator of the SRR oscillation 
amplitudes [Eqs. ^ for ^/-polarization] as well as in all 
elongations in Eqs. ([l2| and Eqs. ([l3|. 

Performing the respective numerical and analytical cal- 
culations as before for the L - structure, we can predict 
the spectral response as well as the effective material 
properties based on the plasmonic eigenmodes. The re- 
sults are shown in Fig. [5] 

Considering these eigenmodes in Fig. [5|d,e), we may 
distinguish two situations. At first an eigenmode which is 
represented by three dipoles being in phase along the en- 
tire structure. The second eigenmode is characterized by 




FIG. 5. (Color online) The far field spectra of the S - structure (a) x-polarization and (b) polarization obtained by numerical 
simulations (circles), predictions based on the SRR structure parameters (dashed-dotted lines), and adapting the parameters 
to fit the numerical values (solid lines). The carrier eigenmodes, i.e., an in line current over the entire structure and antiparallel 
currents [normalized imaginary part of xi(uj) , 1/2(00) , X3(u)] with respect to the center part of the S - structure are observed for 
the two polarization directions (d) x and (e) y. The comparison between the cross-polarization contributions (T^j, Rij) for the 
numerical simulations (circles), the predicted lines from the SRR parameters (dashed-dotted lines) and the fitted parameter 
spectra (solid lines) for the found parameters deduced for fitting the co-polarized response (T^^, Ru). The effective permittivity 
tensor diagonal (c), (f) and the off-diagonal components (i). 



two dipole being in-phase in the arms and out-of-phase 
in the base. Both eigenmodes are excited for x- and y- 
polarization, respectively, and will lead to spectral reso- 
nances appearing as dips and peaks in transmission or re- 
flection, respectively [Fig. |5|a,b)]. The spectral positions 
of the resonances are in agreement with the expectations 
from the SRR structure resonances. 

As expected, for both the S - and the L - metaatom 
the in-phase eigenmodes appear at smaller wavenumbers 
(larger wavelengths) compared to the out-of-phase ones. 
This is completely in agreement with arguments from 
plasmon hybridization theor}^^. 

Again the use of the oscillator parameters as opti- 
mized for the SRR structure [dash-dotted lines in Fig. 
[5|a,b,g,h)] reveals the relationship between both struc- 
tures, since the numerically (circles) calculated spectra 
agree with respect to the overall shape and the resonance 
positions Fig. |5|d,e) very well with the analytical predic- 
tions (solid and dashed lines). A subsequent fitting again 



improves the results towards almost excellent agreement, 
which can also be observed for the cross- polarization 
observables Fig. [5|g,h). Considering the tensor compo- 
nents of the effective permittivity, we observe a differ- 
ence between the diagonal entries due to the geometrical 
differences in the S - structure center and arms, while 
the anti-resonance for the out-of-phase eigenmode is ob- 
served in the off-diagonal elements with the same origin 
as discussed for the L-structure. 

Finally we provide all parameters required for the an- 
alytical calculations presented here for the SRR as well 
as the two presented modifications in Tab. [l| 

V. SUMMARY 

In summary, we have presented an analytical model 
which permits the calculation of effective material pa- 
rameters for planar metamaterials consisting of variable 
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TABLE I. The fitted oscillator parameters applied to repro- 
duce the optical far-field intensities for the SRR as well as 
the adapted parameters for the L- and S- structure are pre- 
sented. For consistence with the corresponding figures, all 
spectral units are given in wavenumber units {cm~^). 



Fitting parameter SRR L S 

cjox,cJoy [cm-'] 9770, 9050 9770, 9500 9770, 9350 

7x,7y [cm-'] 520, 420 520, 420 520, 420 

[cm-'] 6100 6100 6100 

qlr]/eom [AsV/m^Kg] 0.65 • 10^^ 0.83 • 10^^ 0.35 • 10^^ 
qyr]/eom [AsV/m^Kg] 1.1 • 10^^ 0.55 • 10^^ 1.10 • 10^^ 



different metaatoms. Once the constants describing the 
respective dipole oscillations for one particular arrange- 
ment have been determined, where here the SRR has 
been used, the properties of another, modified metaatom 
can be easily predicted. Since in particular the effect of 
asymmetric transmission for circular polarized light at- 
tracted a lot of research interest recently, we focused here 
on planar metaatoms that are optically active. Within 
our model all properties of the effective material tensors 
for such kind of media are correctly predicted and the 
corresponding scattering characteristics are in very good 
agreement with the rigorous numerical results. 



metaatoms formed by a few straight wire sections of po- 
tentially different shape. The model takes advantage of 
resonant electric dipole oscillations in the wires and their 
mutual coupling. The pertinent form of the metaatom 
determines the actual coupling features. Thus this model 
represents a kind of building block approach for quite 
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